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ABSTRACT 

The SuperTIGER (Super Trans-Iron Galactic Element Recorder) instrument was developed to measure the 
abundances of galactic cosmic-ray elements from i 0 Ne to 4 oZr with individual element resolution and the high 
statistics needed to test models of cosmic-ray origins. SuperTIGER also makes exploratory measurements of the 
abundances of elements with 40 < Z ^ 60 and measures the energy spectra of the more abundant elements for Z < 

30 from about 0.8 to 10 GeV/nucleon. This instrument is an enlarged and higher resolution version of the earlier 
TIGER instrument. It was designed to provide the largest geometric acceptance possible and to reach as high an 
altitude as possible, flying on a standard long-duration 1.11 million m 3 balloon. SuperTIGER was launched from 
Williams Field, McMurdo Station, Antarctica, on 2012 December 8, and made about 2.7 revolutions around the 
South Pole in 55 days of flight, returning data on over 50 x 10 6 cosmic-ray nuclei with Z > 10, including ~1300 
with Z > 29 and ~60 with Z > 49. Here, we describe the instrument, the methods of charge identification employed, 
the SuperTIGER balloon flight, and the instrument performance. 
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1. INTRODUCTION 

The SuperTIGER (Super Trans-Iron Galactic Element 
Recorder) instrument, shown in Figure 1 just preceding its flight 
in Antarctica, was developed to measure the elemental compo- 

sition of ultra-heavy galactic cosmic-ray (UHGCR) nuclei with 
charge (Z) 26 ^ Z ^ 40 with high statistical precision and excel- 
lent charge resolution, and to make exploratory measurements 
into the 40 < Z ^ 60 range. A secondary objective was to ac- 
curately measure the energy spectra of the more abundant light 
elements with 12 < Z ^ 30 over the energy range ~0.8-10 GeV / 
nucleon. SuperTIGER was flown during the 2012-2013 Austral 

Summer, returning data on over 50 x 10 6 galactic cosmic-ray 

(GCR) nuclei with Z ^ 10, including ~1300 with Z > 29, and 
~60 with Z > 49 in 55 days at float. 

The origin of GCR nuclei is one of the enduring mysteries 
in astrophysics. Measurements of the elemental and isotopic 
abundances have been used over the last several decades to 
constrain and test theories on the origin of cosmic rays. Casse 
& Goret (1978) and Meyer (1985) suggested that the source 
of cosmic rays was the coronae of stars (types F through 
M). This was based on the observed enrichment of elements 
with low first ionization potential (FIP) in cosmic rays, that 
is similar to that observed in particles from the Sun (solar 
wind and solar energetic particles). According to Schmelz et al. 
(2012), the fractionation in particles from the Sun “probably 
results from a separation of ions and neutrals, which takes 
place between the photosphere and the corona at temperatures 
of 6000-10,000 K.” The coronal material would then be injected 
into the interplanetary medium by stellar winds or flare events 


and later accelerated to cosmic-ray energies by shocks from 
nearby supernovae. In 1989, a paper by Binns et al. (1989, 
p. 997), giving the results of the HEAO-3-C3 experiment 
concluded, “These abundances are consistent with a cosmic- 
ray source having a composition similar to that of the solar 
system, but subject to source fractionation correlated with the 
first ionization potential (FIP) of each element.” 

However, Epstein (1980) and Bibring & Cesarsky (1981) 
showed that refractory elements existing in interstellar dust 
grains should be preferentially accelerated, compared to ele- 
ments existing primarily in the gas phase, owing to their large 
mass-to-charge ratio. The relative enrichment of refractory el- 
ements should be similar to the FIP enrichment since many 
low-FIP elements are refractory and most high-FIP elements are 
volatile. This work was later extended in a pair of papers (Ellison 
et al. 1997 ; Meyer et al. 1997) into a detailed model that reached 
the same conclusion. Additionally, their model predicted that the 
abundances of the volatile elements relative to iron should be 
mass dependent, but refractory abundances should be mass in- 
dependent. Soon after that, Westphal et al. (1998) published 
measurements of the Z = 70-92 charge range from the Trek in- 
strument aboard the Russian space station Mir, which resolved 
even-Z nuclei in this charge range for the first time. They con- 
cluded that the measured relative abundances for Z = 74 through 
82 were inconsistent with a chromospheric origin, but consis- 
tent with an origin in interstellar gas and dust, with refractory 
elements preferentially accelerated. 

More recently, abundances measured by TIGER at 
GeV/nucleon energies (Rauch et al. 2009), ACE at energies 
of hundreds of MeV/nucleon (Binns et al. 2013), and CREAM 
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Figure 1 . SuperTIGER instrument. 

(A color version of this figure is available in the online journal.) 


at TeV/nucleon energies (Ahn et al. 2010), show that refractory 
elements are enriched by a factor of about four over volatiles. 
Furthermore, this enrichment is mass dependent with similar 
slopes for both refractories and volatiles, and the ordering of 
these element abundances with atomic mass is greatly improved 
by comparing them to a mix of normal interstellar material (with 
solar system abundances) and massive star outflow and ejecta 
(Woosley & Heger 2007), rather than normal interstellar mate- 
rial alone (Rauch et al. 2009; Binns et al. 2013). 

The picture that is currently widely accepted is that a large 
fraction of cosmic rays originate in OB associations. The 
OB association origin of cosmic rays was first discussed by 
Montmerle (1979) and Cesarsky & Montmerle (1981). 
More recently Higdon et al. (1998) and Higdon & 
Lingenfelter (2003, 2005, 2013) developed this model in more 
detail. In this picture cosmic-ray nuclei originate primarily 
within OB associations, concentrations of young, short-lived, 
massive stars that form superbubbles in the interstellar medium 
by the combination of their stellar winds and supernova shocks. 
The accelerators are believed to be supernova shocks, shocks 
from the winds of the massive precursor stars (Higdon et al. 
1998; Higdon & Lingenfelter, 2003, 2005; Parizot et al. 2004), 
and colliding winds of massive binaries (Reitberger et al. 2014). 
Cosmic rays are then accelerated from a mixture of old interstel- 
lar gas and dust (with solar-system composition; e.g., Lodders 
2003), the wind outflow from massive stars (including Wolf- 
Ray et stars and their precursor phases), and ejecta from core- 
collapse supernovae (SNII, SNIb,c). The measured composition 
of cosmic-ray isotopes and elements is consistent with a mix of 
interstellar material with outflow and ejecta from massive stars 
(Binns et al. 2005; Rauch et al. 2009), and with refractory ele- 
ments enhanced by about a factor of four over volatile elements. 
Prantzos (2012), on the other hand, can not quantitatively ac- 
count for the measured 22 Ne/ 20 Ne ratio (Binns et al., 2005) in 
his model for an OB association origin. 

Recent observations by the Fermi Large Area Telescope 
(LAT) (Ackermann et al. 2013) of y-rays from two core-collapse 
supernovae, IC443 (a member of the Gem OB-1 association) 
and W44, have confirmed the long-held theory that hadrons are 
accelerated to high energies in at least some supernova remnants 
(SNR). Additionally, distributed high-energy y-ray emission 
has been observed, for example, by Fermi from a “cocoon” 
with morphology coincident with the Cygnus superbubble 
(Ackermann et al. 2011), by VERITAS within the radio-shell 


of SNR G78.2+2 (Aliu et al., 2013), by HESS in the young 
stellar cluster Westerlund 2 (Aharonian et al. 2007), and by 
Milagro in the Cygnus region of the galaxy (Abdo et al. 2007). 
Dermer & Powale (2013) find that 11 of 12 identified y- ray 
emitting SNRs are likely from core-collapse supernovae (SNe). 
Since core-collapse supernovae result from massive stars, and 
most massive stars exist in OB associations, it is clear that a 
large fraction of these SNe occur in OB associations. These 
y- ray measurements are highly complementary to results from 
the direct measurement of cosmic rays and strengthen the OB 
association model of the origin of cosmic rays (Binns 2011). 

TIGER flew over Antarctica in 2001 and 2003 and measured 
the more abundant individual elements up to 33 Sr. The Super- 
TIGER flight significantly improves the statistics of these nu- 
clei, extends statistically significant measurements up to at least 
4oZr, and samples the Z = 50-60 charge range. Sensitive mea- 
surements of these rare Z ^ 30 elements by SuperTIGER will 
enable us to test the OB association/ superbubble origin model 
in more detail by reducing the error bars and extending the reach 
in mass range of the refractory and volatile elements to further 
explore preferential acceleration and its mass dependence for 
ultra-heavy cosmic-ray nuclei. We will also search for possible 
enrichments from nucleosynthesis in massive stars (Woosley & 
Heger 2007). 

As mentioned above, SuperTIGER will also measure the 
energy spectra of the more abundant elements. Heinz and 
Sunyaev (2002) have shown that nearby microquasars may 
impose narrow features in the energy spectra of cosmic-ray 
nuclei in the energy range from ~ 1— 10 GeV/nucleon. The 
complement of sensors used in SuperTIGER allows us to 
measure with greater precision than before (Geier et al. 2006) 
the energy spectra of GCR nuclei up through 30 Zn over the 
energy range of ~0.8-10 GeV/nucleon (top of the atmosphere 
energies). 


2. THE SUPERTIGER INSTRUMENT 

SuperTIGER was designed to maximize the number of 
Z > 30 nuclei that could be measured in a series of long- 
duration balloon flights. This required a design that minimized 
interactions in the instrument and atmosphere by limiting the 
column density of the instrument, using materials with low 
interaction cross-sections, and by flying as high as possible 
for as long as possible. Since long-duration balloon flights 
in Antarctica are limited by the duration of relatively stable 
conditions in the high-altitude winds around the continent, 
SuperTIGER was designed for simple assembly, permitting 
early launches, and used no consumables. SuperTIGER has 
the largest geometric acceptance readily achievable for an 
instrument meeting the weight limits for flight on a standard 
1.11 million cubic meter (40 million cubic foot) balloon. The 
full instrument weight was 1770 kg. 

SuperTIGER incorporates a detector suite selected for ex- 
cellent charge resolution of ultra-heavy galactic cosmic-ray 
(UHGCR) nuclei, minimal nuclear interactions, minimum 
weight, and large geometric acceptance. The techniques used 
to measure charge and energy are identical to those em- 
ployed by TIGER, but the implementation has been improved. 
SuperTIGER (as did TIGER) uses plastic scintillators, acrylic 
Cherenkov counters, and silica-aerogel Cherenkov counters to 
determine the element species and kinetic energy of incident 
particles. Scintillating optical fiber hodoscopes at the top and 
bottom of the detector stack measure the trajectories of parti- 
cles transiting the instrument for angle corrections and mapping. 
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Figure 2. Drawing of SuperTIGER showing the two modules. 
(A color version of this figure is available in the online journal.) 


For SuperTIGER the scintillator and Cherenkov detector enclo- 
sures were floored with an ultra-low-density foam/ aluminum- 
foil composite developed at GSFC, which was also used to 
support the hodoscope (see Section 4 for details). The tops of 
the scintillator, Cherenkov, and hodoscope enclosures were thin 
Al foil. The foam/Al composites and foil windows reduced 
interactions in the structural materials by ~20% compared to 
TIGER. The full geometry factor of both modules combined 
is ~8.3 m 2 sr for particles with trajectory <70° to the zenith. 
Accounting for interaction losses in the instrument, using 34 Se 
nuclei as a reference, the calculated “effective” SuperTIGER ge- 
ometry factor, after interacted particles are removed, is 3.9 m 2 sr, 
~6.5 times larger than TIGER (0.6 m 2 sr). 

To improve detector performance and add redundancy for en- 
hanced reliability, SuperTIGER is divided into two completely 
independent modules, as shown in Figure 2. Each module can 
be mechanically disassembled to a half-module subsystem to 
enable recovery with any available aircraft in Antarctica in- 
cluding the Basler BT-67 and De Havilland DHC-6 Twin Otter. 
Each module has an active area of 1.16 m x 2.4 m, approxi- 
mately the width of TIGER and slightly more than twice as long. 
Each module has its own electronics systems including trigger, 
front-end readout electronics, housekeeping sensors, electronic 
controls, data acquisition/control computer, and data storage. 
SuperTIGER makes efficient use of the active area by mea- 
suring particles that “cross-over” between the active elements 
of the half-modules. SuperTIGER was designed for “graceful 
degradation” and to maintain high performance even with the 
failure of significant numbers of detector channels. 

3. MEASUREMENT TECHNIQUES 

A SuperTIGER module, shown in expanded view in Figure 3 
and fully assembled in Figure 4, consists of three layers of plastic 
scintillator (SI, S2, and S3), a Cherenkov detector with a silica- 
aerogel radiator (CO), another with an acrylic radiator (Cl), and 
two scintillating optical fiber hodoscopes (top and bottom) for 
trajectory determination used to correct for angle of incidence 
and instrument area response. We measure charge and velocity 
with combinations of S and Cl, or Cl and CO, or S and CO. The 
SI and S2 scintillators that make the primary measurements 
of differential energy loss, dE/dx , are located just above the 
top hodoscope (HI) and just below the lower Cherenkov (Cl) 
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Figure 3. Single module shown in expanded view. 

(A color version of this figure is available in the online journal.) 



Figure 4. SuperTIGER module. 

(A color version of this figure is available in the online journal.) 
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respectively. The third scintillation counter (S3) is located below 
the bottom hodoscope (H2), mainly to identify nuclei that have 
fragmented in the instrument, and is also a backup measurement 
in case S2 should fail. Placing H2 between S2 and S3 decouples 
the two scintillators by preventing most 5-rays produced in S2 
from reaching S3. However, we have found that we can obtain 
essentially as good resolution for nuclei with energy below the 
CO threshold by ignoring S3 entirely, by simply demanding 
consistency of SI and S2 signals. 

The use of Cherenkov radiators with different indices of re- 
fraction ( n = 1.043 or 1.025 for CO, and n = 1.49 for Cl) 
enables us to use differing techniques to accurately measure 
charge in complementary energy ranges. For events below the 
threshold energy of CO, charge is measured by the (dE/dx 
vs. Cherenkov) technique using the SI and S2 counters to 
determine dE/dx with velocity corrections from Cl. The use 
of organic scintillators means that the technique is actually 
(dL/dx vs. Cherenkov) where L is the light produced by par- 
ticle energy loss (dE/dx) in the scintillator material. Although 
dE/dx oc Z 2 ,L exhibits saturation effects at high specific dE/dx 
(Birks 1964; Voltz et al 1966; Ahlen et al. 1977; Tarle et al. 
1979; Ahlen 1980; Salamon & Ahlen 1982) and the scintillator 
response becomes more complicated. An approximation for the 
charge dependence of the light output is dL/dx a Z L7 based 
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Acrylic Cherenkov (Cl ) 1 0/20 

Figure 5. Scatter plot of top scintillator signal (S 1 1 / 1 - 7 ) vs. the acrylic Cherenkov 
signal (Cl 1 / 2 ). 
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Aerogel Cherenkov (CO ) 1 0/2 0 

Figure 6. Scatter plot of acrylic Cherenkov signal (Cl 1 / 2 ) vs. the aerogel 
Cherenkov signal (CO 1 / 2 , refractive index 1.04). 


on accelerator calibrations (Binns et al. 1991) and experience 
from TIGER (Link 2003; Rauch et al. 2009). The SuperTIGER 
charge separation in this range is shown in Figure 5, which is 
a scatter plot for a sample of flight data of SI 1/1,7 vs. Cl 0,5 . 
The charge bands are clearly visible. Full correction for scin- 
tillator saturation will use an adaptation of the formulae for 
dL/dx (in Voltz et al 1966; Tarle et al. 1979) based on fits to the 
SuperTIGER data and on a planned calibration at the NASA 
Space Radiation Laboratory at Brookhaven National Labora- 
tory using fragmented 26 ^ Z ^ 54 beams with fragment charge 
determined using precision solid-state detectors. 

For events above the CO threshold, charge is measured by the 
[high-ft-Cherenkov (Cl) vs. low-n-Cherenkov (CO)] technique 
using the Cl signal, with a small velocity correction from the 
CO signal. Above the CO threshold, the Cl signal depends 
weakly on particle velocity. The Cherenkov signals have pure Z 2 
dependence and do not suffer saturation effects. The high energy 
data from SuperTIGER are shown in Figure 6 as a scatter plot of 
Cl 0,5 versus CO 0,5 . Again, the charge bands are clearly visible. 
The [S vs. CO] analysis also provides useful charge separation, 
although resolution is inferior to [Cl vs. CO]. The final charge 
determination is made using a weighted average of the two. 

4. THE INSTRUMENT 

Details of the SuperTIGER instrument and payload systems 
are given below. 


4.1. Scintillation Detectors: 

The scintillators (S-counters) provided both dL/dx measure- 
ments and the event trigger. Each S -counter uses two sheets 
of 1.162 m x 1.162 m x 1 cm ELJEN Technology EJ-208B 
plastic scintillator mounted side by side in a plane (Figure 3). 
This scintillator was chosen for its long attenuation length. The 
scintillation light produced by a nucleus traversing the detector 
is peaked at about 435 nm. This light is transmitted to the edges 
of each sheet by total internal reflection and is coupled through 
a thin air gap to four 1.25 cm square x 1 m EJ-280 wavelength- 


/ Scintillator 
^ Mounting 
Post 



Detector Sides 

PMT Connectors 

Figure 7. Expanded view of scintillator half-module unit. 

(A color version of this figure is available in the online journal.) 


shifter (WLS) bars around its perimeter. The WLS bars shift 
the blue light into the green (peaked at 490 nm); that light is 
transmitted by total-intemal-reflection to Hamamatsu R 1924 A 
2.54 cm diameter photomultiplier tubes (PMTs) coupled to both 
ends of the WLS bars using EJ-500 optical epoxy. The size of 
the WLS bars plus structure result in a 6.4 cm gap between 
the active parts of the side-by-side radiators. A half-detector 
assembly is shown expanded in Figure 7. 

The scintillators were the largest sheets available with the 
required uniformity. The thicknesses of the scintillators were 
measured at the factory using an ultrasonic gauge on an 8 x 
8 grid (15 cm spacing). The flight scintillators were chosen to 
have <0.5%/cm thickness gradient over more than 90% of the 
area with mean gradient <0.25% /cm and maximum gradient 
<1.0%/cm. 

Each scintillator layer within a module was composed of 
two entirely independent light-tight enclosures with custom ex- 
truded Al channel sides and machined corners floored with ultra- 
low density foam/Al composite sandwich. The foam used for 
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the floors was Rohacell™ 31-IG (density 0.032 g cm -1 ), and the 
top and bottom face sheets of the sandwich were 0.1 1 mm thick 
aluminum foil. The thickness of the Rohacell/Al foil sandwich 
was 1.27 cm. The floor of each enclosure is epoxied to the sides 
and corners to form a single, mechanically stable, light-tight 
unit. The enclosure tops are epoxied to 3 mm thick water-cut 
Al frames outside the field-of-view and sealed to the enclosures 
with light-tight custom 1.5 mm thick silicone rubber gaskets. 
The scintillators sit on thin layers of Depron foam, centering 
them on the WLS bars. Sheets of highly reflective front- surface- 
aluminized Mylar (thickness 0.05 mm) above and below the 
scintillators decouple the internal reflection surfaces from the 
foam and improve light collection and signal uniformity. 

The PMT bases used positive high voltage (HV) and radically 
tapered voltage dividers (relative voltage ratios of: K-Dl 1, 
D1-D2 0.17, D2-D3 0.17, D3-D4 0.24, D4-D5 0.24, D5-D6 
0.34, D6-D7 0.34, D7-D8 0.5, D8-D9 0.67, D9-D10 1.0, D10- 
HV 1.0). This gave nearly linear response to signals as large 
as 2 x 10 5 photoelectrons (pe) and an effective dynamic 
range of 2 x 10 4 to cover the charge range 10 ^ Z ^ 
60, including the variations in light reaching each PMT due 
to particle incident position and angle. Typical flight high 
voltage was 650-800 V. The bases incorporated charge-sensitive 
preamplifiers to eliminate pick-up noise that might reduce the 
resolution of the measurement. The PMTs were wrapped with a 
single 0. 1 mm thick layer of mu-metal foil to cancel the effects 
of Earth’s magnetic field which vary with instrument location 
and orientation. 

As in TIGER, the PMTs were mounted using a compliant 
RTV to dissipate vibration during transport as well as to allow 
for thermal expansion/contraction of the WLS bars. However, 
the instrument was exposed to more transport and thermal stress 
than expected before flight, and during integration and test at 
Williams Field in Antarctica a number of the 96 S -counter 
tubes were found to be broken. As a result, we replaced the 
module-2 SI -counter (top counter) with a spare counter. Nine 
additional tubes were not functioning and were presumed to 
also be broken. Because of the difficulty in accessing these 
nine tubes, we elected to cap those tubes off and launch with 
87 functioning tubes. After 18 hr in flight, one additional tube 
stopped functioning. Of these 10 failed tubes, one was in an S-l 
counter, six were in S-2 counters, and three were in S-3 counters. 
These failures did not significantly impact in-flight performance 
due to the SuperTIGER design for graceful degradation. The 
scintillators provided excellent charge resolution below the CO 
threshold and ~100% efficient event trigger for Z ^ 10 nuclei. 

4.2. Cherenkov Detectors: 

The Cherenkov detectors used light integration volumes lined 
with 0.25 mm thick, highly reflective GORE™ DRP light- 
reflective material to collect the Cherenkov light produced by 
the silica-aerogel (CO) and acrylic (Cl) radiators. Each detector 
is a single optical volume 118 cm x 480 cm x 20 cm. For 
ease of recovery it is constructed for easy disassembly into half- 
detectors with dimensions 118 cm x 240 cm x 20 cm. The 
enclosures were constructed of 20 cm x 5 cm Al box extrusions 
machined to accept the PMT housings, with machined end and 
center fittings. The enclosure tops and floors are similar to 
those used in the scintillators. Additionally, the CO aerogel was 
mounted on a Rohacell/Al pallet (Rohacell thickness 0.48 cm 
and 0.008 cm aluminum face sheets) that rested on the floor. A 
full Cherenkov detector is shown with one half expanded and 
the other half assembled in Figure 8. 



Cherenkov Top 


Cherenkov Side Wall 


5" PMTs 
(Hamamatsu R877-100) 


Cherenkov Bottom Panel 


Cherenkov Radiator 
(acrylic or aerogel) 


Figure 8. Cherenkov half-detector showing the assembled (left) and expanded 
(right) module. 

(A color version of this figure is available in the online journal.) 


Each Cherenkov counter is viewed by 42 12.7 cm Hamamatsu 
R877-100 PMTs with “Super-Bialkali” photocathodes for high 
quantum efficiency. These were mounted in custom hydro- 
formed 0.5 mm thick mu-metal sleeves that served as magnetic 
shields, mechanical supports, and light covers. Their bases used 
positive HV with voltage dividers even more radically tapered 
than those for the scintillators (relative voltage ratios of:K-G 1 .0, 
G-Dl 1.0, D1-D2 0.17, D2-D3 0.17, D3-D4 0.33, D4-D5 0.5, 
D5-D6 0.5, D6-D7 0.75, D7-D8 1.0, D8-D9 1.2, D9-D10 1.5, 
D10-HV 1.2). These enabled measurements from 10 pe to ~2 x 
10 5 pe (dynamic range of 2 x 10 4 ) with only 2% nonlinearity 
in order to span the 10 < Z < 60 charge range. As with the 
scintillators, the voltage dividers incorporated integrated charge 
sensitive amplifiers. Flight HV ranged from 950 V to 1200 V. 

Each CO module contained eight aerogel blocks, each approx- 
imately 55 cm x 55 cm x 3 cm. Three of the four half-modules 
contained aerogel blocks with n = 1.043 (12 blocks total), while 
one half-module contained four blocks of index 1.025. These 
have thresholds of 2.5 GeV/nucleon and 3.3 GeV/nucleon, re- 
spectively. The aerogel blocks were purchased from Airglass 
in Sweden in the early 1990’s and were stored in a protective 
dry-nitrogen environment. To improve light production by elim- 
inating any remaining interstitial alcohol and absorbed aerosol 
materials, prior to assembly into SuperTIGER they were baked 
at high temperature in air using a process originally developed 
for IMAX (Mitchell et al. 1996; Labrador et al. 1993; Labrador 
1996). The aerogel blocks are mounted on thin composite pal- 
lets covered in GORE™ DRP and held in position by a layer 
of polyethylene terephthalate, similar to “Saran” Wrap, using 
a technique adapted from BESS/BESS-Polar (Asaoka et al. 
1998). 

Cl uses acrylic radiators with bis-MSB wavelength shifter 
added at a concentration of 25 mg/liter, cast for SuperTIGER 
by Spartech/Polycast. This material has an index of refraction of 
1.49, corresponding to a threshold of 0.3 GeV /nucleon. In each 
module the radiator consists of two 1 16 cm x 1 16 cm x 1 .27 cm 
sheets with a spacing of 4 cm between them. The upper surface 
of the acrylic radiator was blasted with sodium bicarbonate 
powder to produce a light diffusing surface and eliminate total 
internal reflection. The combination of wavelength shifter and 
diffusing surface greatly improves uniformity of light collection 
with position. 

During SuperTIGER development the signals from vari- 
ous enclosure and PMT configurations and radiators were 
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Figure 9. Hodoscope detector showing individual fiber planes (top, x-layer; 
middle, y-layer) and full assembly (bottom). 

(A color version of this figure is available in the online journal.) 


simulated using a Geant4 (http://geant4.web.cern.ch/geant4/) 
Monte Carlo, which uses the measured optical properties of 
the reflector, PMT quantum efficiency, absorption length, and 
Rayleigh scattering length in the Cherenkov radiator. These 
simulations were used to optimize the performance of the light 
collection efficiency, with the result that the charge resolution 
for UHGCR, is superior to that of TIGER even though the Su- 
perTIGER detectors are more than twice as large. The numbers 
of photoelectrons collected corresponded to a mean for singly- 
charged, minimum ionizing particles of 51 photoelectrons for 
Cl, 10 for CO in = 1.043), and 5.2 for CO ( n = 1.025). 

4.3. Hodoscope s : 

Figure 9 is a drawing of the scintillating optical fiber ho- 
doscope. Each module has one x,y - plane on top and one on the 
bottom. Each plane uses two orthogonal layers of square cross- 
section scintillating fibers covering a fiducial area of 2.4 m x 
1.16 m. The short-axis fibers are 1 mm square in cross-section, 
with — 1 160 in each layer. Long-axis fibers are 1 .4 mm square to 
provide improved light transmission efficiency over the 2.4 m 
length, with —828 in each layer. The fibers have polystyrene 
cores and 0.04 mm acrylic cladding. Fluorescent dyes used in 
the core were butyl-PBD (primary dye) at 1% concentration 
by mass and dimethyl-POPOP (secondary dye) at 0.02% con- 
centration. Nuclei with Z = 10 are expected to give >—50 pe 
in a PMT at the greatest distance from the point of incidence 
along the fiber traversed by the particle. This is based on labo- 
ratory measurements of the fibers using a Sr-90 electron source, 
which show that at the farthest distance of the source from the 
PMTs we get >—0.5 p.e.’s. The fibers were drawn over a slowly 
rotating aluminum wheel, then formed into a ribbon glued to- 
gether with Arathane 5753 adhesive and cut to length. The fiber 
fabrication was done in-house at Washington University. The 
ribbons were laid out to form coordinate layers. The fiber layers 
were mounted on a Rohacell/Al composite panel with thickness 
1.59 cm, with the long layer (Y coordinate measurement) glued 


to the panel and the short layers ( X coordinate measurement) 
glued to the long fibers. The short layer is divided into two 
subsections separated by 6. 0-6. 5 cm at the center of the fiber 
layer to approximately correspond to the separation between the 
half-module active areas in the scintillators. 

The hodoscopes use a coded readout scheme originally 
developed for TIGER to limit the number of PMTs required 
(Ward et al. 2011, 2013; Lawrence et al. 1999). The fibers are 
formatted into “tabs” of 6 (long-axis) or 8 (short-axis) fibers, so 
the tab width is ~8 mm. Groups of 12 adjacent tabs at one end 
of each layer are each viewed by a single Hamamatsu R 1924 A 
PMT to give coarse spatial resolution with —9.6 cm segments. 
The tabs on the opposite ends of the fibers are sequentially 
routed to PMTs so each tab within a coarse group is viewed 
by a different PMT. This acts as a “vernier” for single-tab 
localization. A good event has light detected by a PMT at 
both coarse and fine ends. The root-mean- square uncertainty 
in position measurement is o = 8 mm/^/12 = 2.3 mm. With a 
separation of 53 cm between top and bottom hodoscope planes, 
averaged over the effective opening angle of the instrument, 
for Z = 40 the resulting uncertainty in the secant correction 
contributes only cr tra j = 0.05 charge units (cu), which is nearly 
insignificant when added in quadrature with other contributions 
to obtain the overall charge resolution, a z . The PMT signals are 
pulse-height analyzed so that large signals from heavy nuclei 
can be readily distinguished from signals resulting from knock- 
on electrons or accidental protons traversing the fibers. 

Unlike the scintillator and Cherenkov PMTs, a wide linear- 
output range is not required for the hodoscopes. The PMTs used 
conventional linear voltage dividers with flight HV ranging from 
-650V to 950V. 

To account for variations in fiber number and straightness 
between fiber tabs, the X-Y positions of the center of each tab 
on every fiber layer were mapped. Measurements were made 
approximately every 10 cm along the length of the tab to an 
accuracy of 0.5 mm. After the fiber layers were glued to the 
substrate, reference measurements were made to determine the 
absolute position of the fiber layer relative to the substrate and 
frame, thereby locating the position of the center of each tab in 
the instrument coordinate system. These measurements were 
used to determine the positions where tabs intersect on the 
perpendicular layers. These positions were recorded as a lookup 
table allowing the point at which a particle passed through the 
hodoscope to be quickly determined. 

The hodoscope PMT signal for iron nuclei, with trajectory 
<30° to the vertical, as a function of position along the long-axis 
fibers is shown in Figure 10. To determine the attenuation length, 
the x-coordinate was divided into 10 cm bins. Histograms of the 
hodoscope signals from iron events falling into each bin were 
made and fit with a Gaussian distribution. The peak channels 
from the histogram fits are plotted as triangles in Figure 10. 
Fitting a line to these points gives an attenuation length of 
— 1.4 m. The signals falling well below the fitted line are believed 
to be nuclei that traverse significant thicknesses of the acrylic 
cladding, which does not scintillate, instead of the scintillating 
fiber core, or in the case of the drip-lines (e.g., at v-coordinates 
near 200 mm and 650 mm), there may be small gaps between 
the fibers. 

The detection efficiency of the hodoscopes was determined 
by measuring the fraction of events identified as iron-group 
nuclei using the S -counters only for which x,y coordinates were 
obtained at top and bottom. The average detection efficiency 
was 96.5%. 
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Figure 10. Scatter plot of signal in the scintillating fiber long dimension vs. 
distance for iron nuclei. The gap in events near the center of the plot is a result 
of a gap between the two short x-hbers. 

(A color version of this figure is available in the online journal.) 

4.4. Electronics : 

The trigger and readout electronics for a Front End Electron- 
ics (FEE) circuit are shown as a block diagram in Figure 11 
(in this case a Scintillator FEE board). As noted above, each 
of the PMTs in the scintillators and Cherenkov detectors incor- 
porated an integrated charge- sensitive preamplifier. The pream- 
plifiers for the hodoscope detectors were located on the FEE 
boards themselves. The preamp outputs were fed to shaping 
amplifiers with time-to-peak of 1 /xs and read out by peak sen- 
sitive analog-to-digital converters (ADCs). The scintillator and 
Cherenkov ADCs had 16 bit resolution while the hodoscope 
ADCs had 14 bit resolution. For the Cherenkov tubes, two gain 
ranges differing by a factor of 8 were used, each read out by 
16 bit ADCs. 


The scintillator FEE boards and some of the Cl boards 
incorporated fast-shaping amplifiers and comparators. Signals 
registered by the comparators were sent to a field-programmable 
gate array (FPGA) programmed with the trigger logic. Because 
the instrument was only intended to be sensitive to particles 
with Z > 10, the trigger was simple and robust. A logical OR 
was formed of all the PMT signals from each scintillator layer. 
The trigger was logically SI AND (S2 OR S3). Light emitting 
diodes (LEDs) were used to functionally test the PMT systems 
on all detectors. 

4.5. Data Acquisition and Command: 

Each of the two SuperTIGER modules was controlled by its 
own CPU, each with a root disk and a data disk, both based on 
Intel 320 series solid state disks (SSD). The data from both CPUs 
were passed through a multiplexer located on module Ml for 
telemetry through a TDRSS (Tracking and Data Relay Satellite 
System) high-gain antenna that could support transmission up 
to 90 kbps for optimum satellite viewing conditions. A TDRSS 
omnidirectional antenna, with 5 kbps transmission, and an 
iridium antenna were also used. On both modules, an S-band 
transmitter with a downward lobe antenna was used for line-of- 
sight (LOS) communications up to data rates of 150 kbps, when 
the balloon was within line-of- sight range of McMurdo. All 
radio telemetry systems included command uplinks. During the 
flight, initial tuning of the instrument was performed when LOS 
communications were available for both uplink and downlink, 
a period of about 48 hr from the time of launch. After the 
instrument left LOS range, uplink used primarily the Iridium 
link. LOS communications were reestablished for a period at 
the end of the first orbit to increase the downlink rate. 

The average event rate during flight for the combined modules 
was ~30 events s -1 ; the size of each telemetered event (both 
high and low priority) was 398 bytes. The SuperTIGER flight 
software recorded the full flight data on the SSDs, but used a 
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IO BOARD 


Figure 11. Functional block diagram of the scintillator, hodoscope, and Cherenkov readout electronics. Dashed lines indicate additional components of the same kind 
that are not shown. 
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compression scheme for telemetered events whereby only the 
five brightest hodoscope PMTs per side were telemetered (with 
12-bit resolution out of 14 available bits; the remaining 2 bits 
were used to specify the PMT channel number), all Cherenkov 
PMT signals were sent (with 1 bit out of the available 16 used 
to specify whether the PMT signal was in the high or low gain 
channel), while all scintillator PMT signals (16 bit resolution) 
were telemetered. 

All four SSDs failed during flight. Both data disks failed 
at essentially the same time after 10 days at float (Ml on 19 
December at 17:57:47 NZDT and M2 approximately one min 
later, according to the CPU clocks). After that, they would not 
accept write commands. The root disks failed approximately 
10 days later. Fortunately the CPUs did not need to access 
their root disks for data acquisition and telemetry to ground or 
for commanding. Command uplink through iridium and data 
telemetry through TDRSS proceeded successfully for the full 
flight duration. Investigation of the SSD failure has not been 
possible at the time of writing since the instrument has not yet 
been recovered. 

4.6. High-Voltage System : 

To maximize the effective dynamic ranges of the detectors and 
minimize mapping corrections, the high-voltage (HV) system 
was designed to provide great flexibility in tuning the individual 
PMT gains. HV power supplies were connected to groups of 
PMTs through trim boards that allowed each PMT HV to be 
adjusted over a range of ~400 V, providing the ability to balance 
PMT gains. The HV power supplies (HVPS) were compact 
units from EMCO, supplied potted for vacuum use by the 
factory. The Cherenkov counters used both 1250 V CA20P and 
2000 V CA12P supplies, each serving seven PMTs, and CA12P 
were used for the S-counters and hodoscope, each serving four 
scintillator PMTs or six hodoscope PMTs. All connections 
used shielded Reynolds 600 cable. Scintillator and Cherenkov 
PMT cables used Reynolds Series 600 connectors while the 
hodoscope PMT cables were “pig-tails” without connectors. 
The ends of all cables were potted to resist HV breakdown. 
Copper mesh over the HVPS and on the top and bottom of 
the HV trim circuits reduced the local electric field outside the 
potting to help prevent corona. 

The HV systems were tested before flight over a range of 
temperatures (HVPS —22 to +40 C and PMT bases —30 to 
+50 C) and at pressures similar to those expected in flight. 
During flight, one supply had to be turned off due to corona/ 
arcing and five had to be operated at reduced voltage. In addition, 
a number of EMCO voltage monitors showed significant drift, 
but the stable responses of the PMTs indicated that this was 
probably a problem with the voltage monitor signals from 
the supplies, not a problem with the voltages going to the 
PMTs. The failure modes indicated that most of the problems 
were likely due to the EMCO supplies themselves. Although, 
some hodoscope gains were reduced, the instrument remained 
fully efficient for Z > 22 particles and geometric acceptance 
was maintained for the high-priority heavy elements. While 
light collection in one CO detector was reduced, the charge 
determination in the Cl versus CO range depends only weakly 
on CO, and there was no significant loss of performance. 

4. 7. Power System: 

The photovoltaic (PV) power system utilizes two strings 
of 80 SunPower C-60 monocrystalline silicon cells connected 


in parallel. The cells of each string were connected in series 
and laminated onto substrates by SunCat Solar. Each 80-cell 
string consisted of two 5x6 cell arrays and one 4x5 cell 
array mounted in the same plane and connected in series. Each 
contiguous group of 10 cells had a bypass diode installed to 
protect against cell failure and a blocking diode was placed at 
the top of each 80 cell array. The two arrays were mounted in a 
coplanar arrangement on the sun-pointing side of the gondola at 
an angle of about 70° to the zenith (Figures 1 and 2). The other 
PV arrays visible in the photograph on the sides and upper part 
of the gondola, along with those not seen on the other sides of 
the instrument, powered the balloon electronics and telemetry 
systems. These balloon systems were the responsibility of the 
NASA Columbia Scientific Balloon Facility (CSBF). 

The charge controller that the PV arrays fed into was a 
Morningstar TriStar MPPT-45. It was operated in the 24 V 
mode and charged three pairs of Panasonic LC-X1220P lead 
acid batteries. The voltage of the PVs measured at the input to 
the charge controller was very stable at 42 zb 2 V. The system 
provided ~240-260 W of power to the instrument throughout 
the flight. The output of the charge controller was fed into a 
relay circuit board for each module. Relays could be controlled 
by discrete commands. Voltages, currents, and temperatures of 
the power system were monitored. The 24 V main power fed 
DC-DC converters that provided the lower voltages needed by 
the FEE boards, high voltage boards, and CPUs. 

4.8. Gondola, Mechanical Structure, and Thermal Control: 

The size of the instrument components and need to mini- 
mize weight while satisfying CSBF mechanical requirements 
and providing protection for the detectors on landing made the 
gondola especially challenging. Including hodoscope fiber ex- 
tensions and electronics, each module covers an area of 1 .86 m x 
3.1m (Figure 4). Two SuperTIGER detector modules, with spac- 
ing between, required a platform size of 4.3 m x 3.4 m. The 
gondola, shown in Figure 2, uses a light-weight construction 
with commercial Al C-channels forming the deck and the space- 
frame elements made of 3 mm wall Al tubing and machined end 
fittings. 

Rigid 7.5 cm thick honeycomb pallets support each module’s 
detector stacks. These were attached to the gondola deck 
using shock absorbing cable isolators. The detector stacks were 
built up using six C-channel uprights attached to machined 
fittings that were bonded to the honeycomb pallets. This design 
allowed the modules to be fully integrated and tested, including 
electronics, before final integration with the gondola. 

For thermal-vacuum tests at the NASA-Glenn Plum Brook 
facility (2012 June 26-30), the two modules were stacked one 
above the other on a mechanical frame, which also supported the 
modules during shipping to Plum Brook and later to Palestine, 
TX for final integration with CSBF systems. Parts of this frame 
were later used in the shipping boxes that carried the modules 
to Antarctica. 

A comprehensive thermal model was developed for Super- 
TIGER by Scott Cannon at the New Mexico State University 
Physical Sciences Laboratory, and a passive thermal control 
strategy was employed using a single box of foam insulation 
to encase both modules, with 2.54 cm of Techlite insulation on 
top, 5 cm of Styrofoam insulation on the sides and 2.54 cm of 
Styrofoam on the bottom. In flight, most temperatures were in 
the range of —10 to +20° C. However, the top elements of the 
detector, which had low power dissipation, fell to temperatures 
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Figure 12. Flight path of SuperTIGER over Antarctica 2012 December 8 to 
2013 February 1. 

(A color version of this figure is available in the online journal.) 


as low as — 20° C, somewhat lower than the desired minimum of 
— 10°C. The impact of this on the data was minimal. 

5. THE SUPERTIGER FLIGHT 

SuperTIGER was launched on 2012 December 8, from 
Williams Field, Antarctica. In about 2.7 circumnavigations of 
Antarctica, over the trajectory shown in Figure 12, it flew for 
more than 55 days at altitudes from about 36.6 km to 39.6 km 
with an average atmospheric overburden of 4.4 g cm -2 . The 
flight was terminated on 2013 February 1, and set duration 
records for heavy scientific payloads and for heavy-lift scientific 
balloons. The instrument returned excellent data on over 50 x 
10 6 cosmic-ray nuclei above the trigger threshold of ~Z = 10. A 
high-gain TDRSS antenna supported high-rate data downlink, 
up to 90 kbps depending on satellite view angle. A low- 
rate (~5 kbps) TDRSS Omni antenna and an Iridium antenna 
provided command uplinks and backup data links. The latter was 
used for most commanding once the instrument was out of line- 
of-sight (LOS) distance from McMurdo. Good TDRSS links 
were maintained for most of the flight, with typical transmission 
rates of 75 kbps. The SuperTIGER priority system assigned high 
priority to events with signals in the S -counters that exceeded 
that of vertically incident Z = 22 nuclei. The data system 
transmitted high-priority events first and then filled in with lower 
priority data. At 75 kbs, essentially all events were transmitted. 
When poorer links required throttling the transmission to a lower 
data rate, rates as low as 20 kbs still allowed nearly all high- 
priority data to be transmitted. Large fractions of lower priority 
data (10 ^ Z < 22) were also transmitted. Overall, considering 
periods with no telemetry and very-low-rate periods, TDRSS 
data transmission was about 80% efficient for high-priority 
events, with a smaller efficiency for the very abundant low- 
priority events, giving the equivalent of 44 days of data recorded 
for our high-priority data. All triggered events were returned 
during the two LOS (line-of- sight) periods at the beginning 


of the flight and the end of the first orbit/beginning of the 
second. 

The flight was terminated on 2013 February 1 , due to concerns 
about increasing instability in the high-altitude winds and safety 
margins in a CSBF balloon termination battery. Termination 
was carried out “over-the-horizon” and the instrument landed 
on the West Antarctic Ice Sheet (WAIS) at 82.24°S 81.91°W, 
1625 km from McMurdo and ~600 km from the nearest long- 
term US camp (WAIS Divide). Although instrument recovery 
was originally planned for 2014 January, various factors made 
recovery this year impossible, so the recovery is now planned 
for the 2014-2015 austral summer. 

6. PERFORMANCE 

Software that provided housekeeping data displays, scatter 
plots, and PMT histograms, and which incorporated automated 
gain and mapping corrections to the signals, was used to tune 
SuperTIGER after it reached float, by setting HV and thresholds, 
and for monitoring through a web interface during flight. As 
a consequence, the basic performance of the instrument was 
fully verified during the early part of the flight. The methods 
incorporated in this software have been extended to form 
the basis for the instrument analysis software, enhanced by 
additional corrections. 

Following flight, the LOS and TDRSS data were reconciled 
and used to produce a combined data set containing the maxi- 
mum number of “good” events. Data from the hodoscopes were 
then used to determine the straight-line trajectory of each event. 
Pedestal subtraction, PMT gain calibration and adjustments, and 
corrections for detector position-dependent response (mapping), 
were applied. The pedestal and gain values are time and tem- 
perature dependent, and pedestals were measured every 20 min- 
utes. First order charge determinations were made to identify 
iron events, and these were accumulated on a 2 cm x 2 cm 
grid to produce response maps. Corrections were applied to 
each event for PMT time/temperature dependent gain and inci- 
dent angle. The exact positions of the detectors in each module 
were determined with flight data by minimizing the width of the 
Fe peak. 

After initial calibration, the charge (Z) of each particle 
is determined. The procedure is essentially the same for all 
techniques and the [SI vs. Cl] process will be used as an 
example. First, the data are binned by incident angle. Then the 
iron distribution obtained when plotting SI 1/1,7 as a function of 
Cl 0,5 is fit with a quadratic polynomial. (Although dE/dx goes as 
Z 2 , in scintillator the light emitted, dL/dx , goes approximately as 
Z 1 7 due to scintillator saturation; Link 2003). The fit is corrected 
for energy loss through the instrument. For each event, SI 1 / 1 7 
is scaled by the corrected fit, placing the events on a charge 
scale with iron at Z = 26. The assigned-charge data for all 
angle bins are combined and residuals in the iron curve as a 
function of Cl are fit with a quartic polynomial. The charges 
are corrected using this fit to produce a final charge assignment. 
All assignments using the scintillators employ a corresponding 
procedure with fits as a function of either Cl 0,5 or CO 0,5 . In the 
[Cl vs. CO] technique, the fits are made in Cl 0,5 as a function 
of CO 0,5 . 

The data are analyzed as two separate data sets; events that 
have a CO signal above threshold and those that do not have 
a CO signal above threshold. Below the CO threshold, charges 
are determined separately using [SI vs. Cl] and [S2 vs. Cl] 
as described above. The charge assigned is a weighted average 
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Charge 

Figure 13. Histogram of events in the energy range below the aerogel Cherenkov 
threshold (2.5 and 3.3 GeV/nuc for refractive indices ( n ) = 1.04 and 1.025, 
respectively). The charge was determined using SI, S2, and Cl. 

(A color version of this figure is available in the online journal.) 


of the two. The resulting charge histogram for 10 ^ Z < 31 
is given in Figure 13. For these data, agreement was required 
between SI, S2, and S3, but only SI, S2, and Cl were used 
to determine charge. The Gaussian fit to the Fe peak gives 
a resolution in charge of oz = 0.16 cu, demonstrating the 
excellent performance of SuperTIGER. This compares with the 
oz = 0.23 reported by TIGER (Link et al. 2003; Rauch et al. 
2009). Note that odd-Z 27 C 0 and 29 CU are clearly visible in the 
histogram even though they are adjacent to the ~ 100 times more 
abundant 26 ^e and 28 Ni peaks. Note also the clearly resolved 
3 oZn peak. 

Above the CO threshold, charge is determined using a 
weighted average of [Cl vs. CO] and [SI vs. CO]. As was the 
case for the events with CO signal below threshold, agreement 
was required between SI, S2, and S3 signals. The resulting 
charge histogram for 10 ^ Z < 31 is given in Figure 14. The 
Gaussian fit to the Fe peak has oz — 0.16 cu, as was the case for 
the “below-CO” events. This data set also shows clearly visible 
27 Co and 29 Cu peaks and a clearly resolved 3 oZn peak. 

The number of iron events in these above- and below-CO data 
sets is 2.98 x 10 6 . However, the data set can be expanded to 
obtain a larger statistical sample for the very rare ultra-heavy 
nuclei by relaxing the agreement criteria. In this expanded data 
set (not shown), for the case of events below the CO threshold, 
we require agreement between SI and S2 only (S3 is not used 
in the analysis at all), and for events above the CO threshold 
we require agreement only between SI and Cl (S2 and S3 are 
not used in the analysis). This data set contains 3.98 x 10 6 
iron nuclei, with nearly the same charge resolution at iron of 
about 0.17 cu. For comparison the number of iron events in 
the combined 2001 and 2003 TIGER flights was 0.56 x 10 6 
(Rauch 2008), giving a factor of 7.1 increase in the number 
of analyzable events over that of TIGER. The data set also 
contains 1130 nuclei with Z > 30.5 and ~60 nuclei with 
Z> 49.5. 

At the present stage of analysis clear peaks can be observed up 
to 44 Ru (not shown). Work is ongoing to improve the resolution 
in this charge range, and to extend it up through the Z = 
50-60 range (with a low statistical sample) through the use 
of more sophisticated charge and velocity dependent scintillator 
saturation corrections as discussed in Section 3. 



Charge 

Figure 14. Histogram of events in the energy range above the aerogel Cherenkov 
threshold (2.5 and 3.3 GeV/nuc for refractive indices (n) = 1.04 and 1.025 
respectively). The charge was determined using Cl, CO, and SI. 

(A color version of this figure is available in the online journal.) 

7. CONCLUSION 

The SuperTIGER instrument was developed to accurately 
measure the abundances of ultra-heavy, trans-iron GCR nuclei. 
It performed well during its record 5 5 -day Antarctic flight, 
returning data on over 50 million cosmic-ray nuclei. Initial 
analysis of these data has shown that the charge resolution of 
the instrument for iron nuclei is 0.16-0.17 charge units at all 
energies. The numbers of iron nuclei that we have detected 
is ~4 x 10 6 , a factor of 7.1 larger than for the combined 
TIGER flights in 2001 and 2003. Work is ongoing to extend 
this excellent resolution into the Z = 50-60 range. 
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